On the nature of spinor Bose-Einstein condensates in rubidium 
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We perform detailed close-coupling calculations for the rubidium isotopes S5 Rb and 87 Rb to 
ascertain the nature of their spinor Bose-Einstein condensates. These calculations predict that the 
spinor condensate for the spin-1 boson 87 Rb has a ferromagnetic nature. The spinor condensates for 
the spin-2 bosons 85 Rb and 87 Rb, however, are both predicted to be polar. The nature of a spin-1 
condensate hinges critically on the sign of the difference between the s-wave scattering lengths for 
total spin and 2 while the nature of a spin-2 condensate depends on the values of the differences 
between s-wave scattering lengths for the total spin 0, 2 and 4. These scattering lengths were 
extracted previously and found to have overlapping uncertainties for all three cases, thus leaving 
the nature of the spinor condensates ambiguous. The present study exploits a refined uncertainty 
analysis of the scattering lengths based on recently improved result from experimental work by 
Roberts et al, which permits us to extract an unambiguous result for the nature of the ground state 
spinor condensates. 
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In a conventional magnetic trap for ultra-cold alkali 
atoms the spin degrees of freedom are "frozen out" since 
the atom must be in a weak-field seeking Zeeman state to 
be trapped. In an optical trap, however, the spins of the 
alkali atoms are essentially free, and all magnetic sub- 
states |/, to) for a given spin / can be populated. Since 
the atom-atom interaction depends on spin, these mag- 
netic substates can be changed in a scattering event. Ac- 
cordingly, it is of interest to see how the spins are orga- 
nized in the ground state and to explore the nature of 
the spin-mixing dynamics in an optically trapped Bose- 
Einstein condensate. 

Multi-component condensates have been formed in 
magnetic traps. For instance Ref. |l| used a dou- 
ble magneto-optical trap and a magnetic trap to cre- 
ate condensates in cither the \f — 2, to = 2) or the 
|/ = 1,to = —1) spin state of 87 Rb, and in a mixture 
of both by cooling |1, — 1) evaporatively and |2,2) via 
thermal contact with the |1, — 1) atoms. In this case the 
spin projections are approximately frozen out because 
the spin-flip cross sections in 87 Rb are anomalously small 
|^|, ||, By contrast, Ref. Q made a sodium condensate 
consisting simultaneously of all three magnetic substates 
of the / = 1 atomic state, by cooling the atoms in a 
magnetic trap and then transferring them into a optical 
trap. This experimental technique produces what is re- 
ferred to as a spinor condensate, because it can explore 
its full range of spin degrees of freedom. In the theo- 
retical description of Ref. || M, the spinor condensates 
are classified according to the relative values of certain 
characteristic scattering lengths. Note that alternative 
theoretical treatments |j| differ in their detailed pre- 
dictions concerning the nature of the spinor BEC ground 
state. Nevertheless, in this paper we determine the in- 
teraction parameters for spinor condensates of 85 Rb and 
87 Rb which based on Ref. §, @, fall into the two 
following categories: 

Spin-1 atoms ( 87 Rb) Let F be the total spin of two 



bosonic spin / = 1 atoms, and let clf be the s-wave 
scattering length for the total spin F symmetry. Since 
/ = 1, only F = 0, 2 are allowed by Bose symmetry 
for an s-wave collision. The nature of the spin-1 BEC 
ground state depends critically on the relative values of 
ao and 02- According to Ho || a spinor Bose conden- 
sate composed of spin-1 bosons in an optical trap can 
be either "ferromagnetic" or "antiferromagnetic" in na- 
ture H 0] . The antiferromagnetic state has alternatively 
been termed "polar", and we use this terminology here. 
The difference between the scattering lengths ao and ai 
determines the nature of the spin-1 condensate: the ferro- 
magnetic state emerges when ao > 0,%, whereas the polar 
state emerges when ao < ai In the ferromagnetic 
state virtually all atoms reside in the same spin substate 
(either to = 1 or m = — 1); in the polar state the spin 
projections are mixed. 

Spin-2 atoms ( 85 Rb, 87 Rb) Two bosonic spin / = 2 
atoms possess F — 0, 2, 4 total spin states exhibiting 
the appropriate Bose symmetry for an s-wave collision. 
For spin-2 87 Rb the scattering lengths ao, a 2 and 04 are 
determined by the real part of the phase-shift since the 
inelastic scattering processes are also allowed. According 
to Ciobanu et al. pToj ] , a spinor condensate of spin-2 
bosons in an optical trap can be one of the three types 
"ferromagnetic", "polar", or "cyclic" in nature, which 
we abbreviate as F, P or C respectively. Ferromagnetic 
and polar condensates are similar to those above. The 
name "cyclic" arises from a close analogy with d-wave 
BCS superfluids. The nature of the spin-2 BEC ground 
state depends critically on the relative values of ao — 02 
and a 2 — 04 Jl(| . 
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The three states emerge under following conditions: 



P: a - a 4 < 0, f \ a 2 - a 4 \< ± | a - a 4 |, 
F: a 2 - a 4 > 0, \ | a - a 4 | +|(a 2 - a 4 ) > 0, 
C: a 2 — a 4 < 0, -| | ao — a 4 | — |(a 2 — a 4 ) > 0. 
For spin-1 Rb the total spin F = and F — 2 scat- 
tering lengths do and a 2 are almost equal. They have 
been calculated before |J based on the analysis of Ref. 
jllj , but the uncertainties determined still overlap for ao 
and a 2 , so that the sign of the difference has remained 
uncertain. In particular, the scattering lengths have been 
interpreted rather conservatively in Ref. [[flj. For spin- 
2 85 Rb and 87 Rb the uncertainties for the total spin F 
scattering lengths ao, a 2 and a 4 have been too large to 
uniquely identify the nature of the spinor condensates 
M. The uncertainty region for 85 Rb was large enough to 
overlap all three regions P, F and C, while the uncertainty 
region for 87 Rb overlapped both the polar and the cyclic 
region. In the present study we determine the scatter- 
ing lengths ao and a 2 , and their uncertainties for spin-1 
87 Rb, and ao, a 2 and a 4 and their uncertainties for spin-2 
85 Rb and 87 Rb. We concentrate on an accurate determi- 
nation of the difference ao — a 2 for spin-1 87 Rb and the 
pair (ao — a 4 , a 2 — a 4 ) for spin-2 85 Rb and 87 Rb. If one 
accepts the spinor condensate treatment of Ref. || ||] 
this analysis gives an unambiguous determination of the 
nature of the BEC ground states. 

Uncertainties in the scattering lengths arise primarily 
from imperfect knowledge of three parameters: the long- 
range van der Waals coefficient Cq, and the singlet and 
the triplet s-wave scattering lengths, a s and a t respec- 
tively. In addition, when using potential curves deter- 
mined for one isotope to predict scattering for another 
isotope, the results can depend on the precise number of 
bound states in the triplet potential, Nb as well as the 
precise number of bound states in the singlet potential. 
Roberts et al. analyzed a magnetic-field Feshbach reso- 
nance to determine "state of the art" potentials for 85 Rb 
p| . Recently they have revisited some of the rether- 
malization measurements in Ref. 11 and improved the 
uncertainties for the long-range van der Waals coefficient 
and the singlet end triplet s-wave scattering lengths for 
85 Rb |l|. 

Using these new values of Cq, a s and at we show below 
unambiguously that ao > a 2 for spin-1 87 Rb. This result 
in turn implies that the spinor condensate is definitely 
ferromagnetic, as was previously suspected ||. By con- 
trast the spin-1 23 Na scattering lengths, recently deter- 
mined in Ref. Jl3| ], imply that a 23 Na / = 1 spinor BEC 
is polar, as has been suggested before ||. By extracting 
the scattering length from a spectroscopic experiment, 
Crubellier et al. found that, for 23 Na, ao = 50.0±1.6 a.u. 
and a 2 = 55.0 ± 1.7 a.u. They calculated the scat- 
tering lengths for two values of the Cq coefficient for 23 Na 
and found that the influence of the Cq value is very small 
(a 4% change in Cq results in a variation in the scattering 
length of the order of 0.1%). Consequently, the analysis 
for 23 Na Jl3t , in conjunction with the present analysis for 
87 Rb, implies that both types of spin-1 condensates can 



be realized with the atoms used most frequently in BEC 
experiments ( 23 Na and 87 Rb). 

The improved results for Cq , a s and a t also predict that 
ao — a 4 < 0, | | a 2 — a 4 |< -i | ao — a 4 |, for both spin-2 
85 Rb and 87 Rb. This result implies that both spin-2 85 Rb 
and 87 Rb will be polar. Previously, it was estimated that 
87 Rb would be polar, but that 87 Rb would be cyclic ]lO| |. 
This implies that the ground state for spin-2 85 Rb and 
87 Rb will have the same nature as spin-2 23 Na. Spin-2 
23 Na was already unambiguously classified since the un- 
certainties on differences between the relevant scattering 
leng ths place ao — a 4 and a 2 — a 4 within the polar region 
Jlbf, The results for spin-1 87 Rb and for spin-2 85 Rb and 
87 Rb are summarized in Figs. |l], ^| and |3| respectively. 

Our calculations start from the singlet and triplet 
Born-Oppenheimer potentials between two rubidium 
atoms that were calculated in Ref. |Q, where the sing- 
let potential is adjusted to have 125 bound states p5[ . 
These potentials are matched smoothly at r = 22.0 a.u. 
to the standard long-range van der Waals potentials us- 
ing the new value of the long-range coefficient Cq in- 
ferred from the experiment in Ref. |Tl| and reanalyzed 
according to Ref. p2| , and using the C 8 and C w coeffi- 
cients from the calculations of Ref. jl(| . The potentials 
are adjusted to match the scattering length by includ- 
ing short-range inner-wall corrections that are parame- 
terized for each spin by: carctan((r — r m i n ) 2 / (c,-)) for 
T < r m i n ■ c r is a constant (the same order of magnitude 
as Train, slightly different for the singlet and the triplet), 
the inner-wall parameters c are of the order of 10~ 5 to 
10~ 4 a.u., r is the separation between the two Rb atoms 
and r m in is the separation for which the potential is min- 
imal. The inner- wall parameters c are varied over a range 
that reproduces the recently improved values of a s and 
a t . The improved values of Cq for rubidium, a s and a t 
for 85 Rb are: Cq = 4660 ± 20 a.u., a s = 3650le7o a.u. 
and at = —332 ± 18 a.u. |Q, while the calculations of 
Ref. |T(| determined that Cg = 550600 a.u.. These are 
the values we adopt in the present calculations. Our cal- 
culations here do not allow for variance in C%- This is 
reasonable because the dependence of Cg is one order of 
magnitude smaller than the dependence of Cq. Further- 
more the number of bound states in the triplet potential 
was previously believed to be 39 ± 1 (lj, [TtJ , but more 
refined experimental analysis suggests that it is instead 
40 < N b < 42 JT|, [l8| . The present calculations are done 
for N b = 39,40,41,42. The number of bound states in 
the singlet potential is not changed. 

Our calculations have been carried out for three values 
of Cq that span the empirical range (4640 — 4680 a.u.). 
These values are adequate since the quantities of interest 
vary smoothly with Cq over the range of interest. We 
also tested the triplet potential for each one of the four 
relevant N b (39, 40, 41, 42). For each value of N b we 
determine the values of the inner-wall corrections that 
correspond to the uncertainty range of 85 Rb a s and at 
for each of the three values of Cq. These calculations are 
carried out at zero magnetic field and 130 nK since the 
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FIG. 1: Total spin F = scattering length versus total spin 
F = 2 scattering length for spin-1 87 Rb. The uncertainties 
of ao and a 2 are determined by the uncertainties on a s , at, 
Ce and Nb, the number of bound states in the 85 Rb triplet 
potential. The symbols in the middle of the "diamonds" are 
the mean scattering length for each Ce and the diamonds 
encircle the uncertainties arising from uncertainties on a s and 
at for all Ce- The thick black line shows the boundary between 
the ferromagnetic and polar phase of the spinor condensate: 
ao = a 2 and the number next to each diamond is Nf,. 



given value of Ce, a s and a t are determined from colli- 
sions at this temperature [T^| . The same potentials opti- 
mized for 85 Rb have been used in our 87 Rb calculations, 
except for an appropriate change in the reduced mass. 
Since Nb in rubidium is unknown at present, and since 
we utilize the same potentials determined by the 85 Rb 
singlet and triplet scattering length in our 87 Rb calcula- 
tions, we incorporate N b in our analysis of the uncertain- 
ties for 87 Rb. The singlet and triplet potentials are used 
in multichannel calculations to compute the ao and 02 for 
spin-1 87 Rb and ao, a2 and 04 for spin-2 85 Rb and 87 Rb 
(again at zero magnetic field and at 130 nK), with fixed 
Ce and Nb- The calculations are then repeated for each 
value of Cq and Nb with the corresponding new values of 
the inner- wall corrections. These calculations span the 
empirical Ce, Nb, a s and a* range, which permits us to 
extract the over-all uncertainty in the difference between 
the relevant scattering lengths, ao — a2 for spin-1 and the 
two relevant differences (ao — 04, ai — a.4) for spin-2. 

As a confirmation, the same rubidium potentials have 
also been used to calculate the singlet and triplet s-wave 
scattering lengths a s and at for 87 Rb (at 130 nK energy). 
These single-channel calculations have been repeated for 
the three values of Ce and four values of Ce that span 
the empirical range. This permits us to check whether a s 
and a t fall within the range of previous measured values. 
The single-channel triplet scattering lengths are found to 
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FIG. 2: Difference between total spin F = 2 and F = 4 scat- 
tering lengths versus the difference between total spin F — 
and F = 4 scattering lengths for spin-2 87 Rb. The uncertain- 
ties of 00, 02 and a^ are determined by the uncertainties of a a , 
at, Ce and Nb, the number of bound states in the 85 Rb triplet 
potential. The symbols in the middle of the "diamonds" are 
the mean scattering length for each Ce and the diamonds en- 
circle the uncertainties arising from uncertainties on a s and 
at for all Ce- The thick black line shows the boundary be- 
tween the polar and the cyclic phase of the spinor condensate: 
(oi2 — 0.4) = ^(a2 — an) and the number next to each diamond 
is Nb- 



be: 

N b a t 

39 107 ±1 a.u. 

40 103 ±1 a.u. 

41 100 ±1 a.u. 

42 97 ±1 a.u. 

and the single-channel singlet scattering length is found 
to be: a s — 91 ± 1 a.u.. The a s and at values for Nb = 39 
are in good agreement with previous work ]TT] , |i"9| . As 
another confirmation, and since Nb is unknown we have 
also calculated the 85 Rb scattering length for / = 2, m = 
—2 at various magnetic fields to compare the obtained 
values with the values from Roberts et at; ]l2[ ]. This 
comparison shows good agreement. For each of the A^, 
the scattering lengths obtained for the specific magnetic 
fields exhibits an uncertainty greater than the one given 

by @. 

The values for ao and 02 for spin-1 87 Rb, along with 
their uncertainties, are shown in Fig. [l]. ao is always 
greater than a-i in the multichannel calculations, which 
unambiguously determine the nature of spin-1 87 Rb to 
be ferromagnetic. The global difference lies between 0.3 
and 2.7 a.u. over the uncertainty range. The difference is 
an increasing function of Nb, while ao and a 2 themselves 
are decreasing functions of Nb. For a given Nb, the range 
of possible values of ao and a2 varies only weakly with 
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FIG. 3: Difference between total spin F = 2 and F = 4 scat- 
tering lengths versus the difference between total spin F — 
and F — 4 scattering lengths for spin-2 85 Rb. The uncertain- 
ties of arj, 02 and 04 are determined by the uncertainties of a s , 
at, C% and Nb, the number of bound states in the 85 Rb triplet 
potential. The symbols in the middle of the "diamonds" are 
the mean scattering length for Nt, = 41 and each Ce and the 
diamonds encircle the uncertainties arising from uncertain- 
ties on a s and at for each Ce- The thick black line shows 
the boundary between the polar and the cyclic phase of the 
spinor condensate: (02 — 04) = -^((12 — a-i) and the close-up 
shows the average values for the four different Ni,. 

C$, as was the case for 23 Na |l3| . 

The results of an — a 4 and a 2 — 0,4 with uncertainties for 
87 Rb spin-2 are shown in Fig. |[ For all four values of Nb 
and all three values of Cg the uncertainty region is within 
the "polar" region, making the nature of 87 Rb spin-2 
condensate unambiguously determined. The pair (an — 
04, a 2 — 04) moves closer to the boundary between the 
polar and the cyclic regions as Nb is increased but never 
reaches the boundary, within the present uncertainties. 
For a fixed value of Nb ai — 04 is increasing as a function 
of Cg, while an — 04 is almost independent of Cg. The 
uncertainty region for a fixed value of Nf, is very narrow 
(especially for higher Nb). The long axis of this region 
corresponds to the difference a s — at, whereas the narrow 
axis corresponds to the sum a s + a t . 

The results for 85 Rb spin-2 are shown on Fig. § In 
contrast to the case of 87 Rb, here a-i — 04 and an — 04 are 
more dependent on the value of Cg than on Nf,, but only 
very little. 02 — 04 and an — 04 are slowly increasing func- 
tions of Cg as well as of Nb. The uncertainties for all Nb 
and values of Cg unambiguously determine the nature of 
85 Rb spin-2 to be polar. The uncertainty region is again 
very narrow. The long axis of this region corresponds to 
at, whereas the narrow axis corresponds to a s . 

Since the graphs for spin-2 85 Rb and 87 Rb show 



only scattering lengths differences rather than scatter- 
ing lengths, we summarize an, a2 and 04 in the following 
table. The scattering lengths for 85 Rb show only very 
little dependence of Nb. Over the entire range of a s , a t , 
Cg and Nb the estimated scattering lengths for 85 Rb are 
(in a.u.): 

a Q = -680 ± 80, a 2 = -540 ± 50, a 4 = -390 ± 30 
The scattering lengths for spin-2 87 Rb over the entire 
range of a s , a t and Cg are estimated as (in a.u.): 

N b a a 2 04 

39 90.5 ±1 97.5 ±1 106.8 ± 1 

40 89.0 ±1 94.8 ±1 103.6 ±1 

41 87.7 ±1 92.4 ±1 100.5 ± 1 

42 86.4 ± 1 90.2 ± 1 97.4 ± 1 

These numbers conservatively give the global uncertain- 
ties for each Nb. In the context of spinor condensates it 
is necessary to consider the actual allowed regions of the 
parameters, as we have done above and which permit us 
to draw meaningful conclusions. 

To see how the results change when the multichan- 
nel and single-channel energy changes, we calculated the 
scattering lengths at various energies (with 1 pico Kelvin 
in the single and multichannel as the lowest value) to 
cover the relevant temperature for some experiments. 
This did not change our conclusions about the nature 
of the spinor Bose-Einstein condensates in rubidium. 

Since the shape of the inner-wall potential is not known 
exactly and since we change it to have potentials with the 
four different values of Nb we also performed the calcula- 
tions with a quadratic inner- wall correction (c(r — r m i n ) 2 
for r < r m in) instead of the arctan form. This did not 
change the conclusions and only changed the calculated 
scattering lengths by about 0.1 %. 

The present values of an, ai for spin-1 87 Rb and an, 
a%, 04 for spin-2 85 Rb and 87 Rb are consistent with val- 
ues obtained from Ref. jnj. Note that to carry out the 
calculations based on Ref. jTTJ , the correlations among 
a s , a t and Cg must be taken into account. We have sep- 
arately calculated the values of an, 02 for spin-1 87 Rb 
and an, ai, 04 for spin-2 85 Rb and 87 Rb form a s , at for 
85 Rb and Cq as given in Ref. |ll[], and find that they 
support our classifications of the spinor condensates as 
presented in this paper. The new values from Ref. |l2| 
allow us to determine a smaller uncertainty on the cal- 
culated scattering lengths, but they do not change our 
conclusions. 

In summary our analysis based on the new results for 
the values of Cq, a s , at and the number of bound states 
in the triplet potential demonstrate that the nature of 
the ground states of 85 Rb and 87 Rb spin-2 condensates 
should be polar. In addition, the ground state of the 
87 Rb spin-1 condensate should be ferromagnetic. There- 
fore, in view of the known scattering parameters for 23 Na, 
both ferromagnetic and polar spin-1 condensates are ex- 
perimental accessible whereas no cyclic or ferromagnetic 
spin-2 condensate appears to exist for the most common 
rubidium isotopes. 
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